The spin-triplet superconductivity induced by the charge fluctuation in extended 

Hubbard model 
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The pairing symmetry in the electron mechanism for superconductivity is explored when charge 
fluctuations coexist with spin fluctuations. The extended Hubbard model is adopted to obtain, with 
the fluctuation exchange approximation, a phase diagram against the on-site Coulomb repulsion 
U and the off-site repulsion V for the square lattice with second-neighbor hopping t . We have 
found that (i) for large U(> 9) a triplet superconductivity with a sin(k x + k y ) symmetry can appear 
just below the charge density wave phase. The pairing is degenerate with sin(k x — k y ), so a chiral 
sin(k x +k y )-\-i sm(k x —k y ) that breaks the time reversal symmetry should result, which is a candidate 
for the gap function on the 7 band of Sr2Ru04 and is consistent with a recent measurement of the 
specific heat, (ii) By systematically deforming the Fermi surface with varied t' , we have identified 
the region where the triplet pairing is most favored to be the region where the Fermi surface traverses 
the van Hove singularity with the charge susceptibility strongly enhanced. 

PACS numbers: 74.20.Mn 



I. INTRODUCTION 



Spin-triplet superconductivity, which is arousing much 
interests in recent years, is fascinating in a number of 
ways. Theoretically, a most intriguing question is the 
following: if we consider the electron mechanism of su- 
perconductivity in the most frequently adopted Hubbard 
model with an on-site electron-electron repulsion, we can 
show that triplet superconductivity is very difficult to 
realize for the simple reason that the pairing interaction 
mediated by spin fluctuations is only 1/3 in the triplet 
channel than in the singlet channel 1 ' 2 . So any theory at- 
tempting to explain triplet pairing has to overcome this 
question. 

Experimentally, the discovery of superconductivity by 
Maeno and coworkers in the layered perovskite ruthe- 
nium oxide Sr2RuC>4 have kicked off renewed inter- 
ests. Suggestions for a triplet pairing in this mate- 
rial came from NMR Knight shift 3,4 and polarized neu- 
tron scattering 5 . A broken time reversal symmetry is 
further observed with /xSR 6 ' 7 and small-angle neutron 
scattering 8 . As for the pairing symmetry, NMR and 
NQR relaxation rates have shown an absence of the 
Hebel-Slichter peak 9,10 , which suggests line node(s) in 
the gap function, which was supported by a specific heat 
measurement 11 . While the angular dependence of ther- 
mal conductivity indicates the presence of horizontal line- 
nodes 12 , a recent field-orientation dependence of the spe- 
cific heat shows that the gap in the active 7 band has 
minima along [100] directions with the passive a and (3 
bands having gap minima along [110] directions 13 . 

A theoretical work by Zhitomirsky and Rice 14 has pro- 
posed horizontal line nodes with a good fit to the specific 
heat measurement 11 . As for the mechanism that sta- 
bilizes triplet pairing, Kuwabara and Ogata 15 , and in- 



dependently Sato and Kohmoto 16 , have suggested that 
triplet p-wave pairing can be induced by an anisotropy 
in the spin fluctuation. Kuwabara and Ogata 15 have 
identified the most suitable gap function to be sin k x + 
i sin k y that breaks the time-reversal symmetry, which 
was also proposed by Miyake and Narikiyo 17 . The chiral 
<Ar(fc) + { < t>x(k),4>y(k): p-wave like) has also been 

obtained with a third-order perturbation theory by No- 
mura and Yamada 18,19 and more recently by Yanase and 
Ogata 20 , and with the fluctuation exchange approxima- 
tion (FLEX) for quasi-one-dimensional a and (3 bands by 
Kuroki et al. 21 . 

However, the validity of perturbation results truncated 
at finite orders has to be checked by non-perturbative 
methods. Indeed, a recent quantum Monte Carlo (QMC) 
result by Kuroki et al. have shown that the singlet im- 
pairing dominates over the triplet for the 7 band. 22 This 
situation has led Arita et al. to show, with the dy- 
namical cluster approximation (DC A), that d x 2_ y 2 pair- 
ing actually dominates over p-wave pairings as far as 
the on-site Hubbard model is concerned, while if we go 
over to the extended Hubbard model with off-site repul- 
sion the triplet superconductivity can be favored but not 
dominant in extended Hubbard model employing nearest 
neighbor Coulomb repulsion for 7 band in Sr 2 Ru04. 23 

Generally, there is a FLEX result by Arita et al 1 who 
show that triplet superconductivity is much weaker than 
singlet pairs for the one-band, on-site Hubbard model, 
which agrees with a phenomenology 2 , while if charge 
fluctuations are enhanced by off-site repulsions triplet 
pairs have a chance to dominate. We have already 
shown triplet superconductivity can be dominant near 
the charge density wave (CDW) phase on square lattice 
with the FLEX in the extended Hubbard model. 24 

So the purpose of this paper is to examine how triplet 
superconductivity can become dominant in the extended 
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Hubbard model with the FLEX in general, and for 7 
band of Sr2Ru04 in particular. The extended Hub- 
bard model has been studied, primarily for specific 
charge densities n, e.g., half-filling or quarter filling, by 
means of the quantum Monte Carlo method 25 , the weak 
coupling theory 26 , the mean-field approximation 27 ' 28 , 
the second-order perturbation 29 , the random phase 
approximation 30 , the FLEX approximation 31 , the slave- 
boson technique 32 , the bosonization and the renormaliza- 
tion group 33 ' 34 . Here we adopt the FLEX developed by 
Bickers et a/. 35-38 , which is a renormalized perturbation 
method to study pairing instabilities when exchange of 
spin and charge fluctuations are considered as dominant 
diagrams. This approximation is useful to explore the 
tendencies of dominant pairing relatively when the pa- 
rameters are systematically varied. However, including 
the off-site Coulomb repulsion V susceptibilities and ef- 
fective interactions become (Z + l) x (Z + l) matrices for 
the lattice coordination number Z (= 4 for the square lat- 
tice), which demand extremely computer resources. We 
show as far as the present finite-temperature result is 
concerned that (i) triplet sin(k x + k y ) pairing appears in 
between singlet cos 2k x — cos 2k y and the CDW phase for 
U > 9, and (ii) the pairing symmetry changes as d x i_ y i 
— > sm(k x + k y ) — > cos 2k x — cos 2k y when the shape of the 
Fermi surface is varied by the second-neighbor hopping 
t'. Physically, all the results can be well explained by the 
structure and peak value of spin and charge susceptibili- 
ties. 



II. FORMULATION 

Let us start with the extended Hubbard Hamiltonian, 

nn.nnn 



tin 



ij TlifjTljtj' , 



i.j crcr' 



in the standard notation on a tetragonal lattice depicted 
in Fig.l, where nn (nnn) denotes nearest neighbor (next- 
nearest neighbor) sites. For the square lattice the unit of 
energy is taken to be the nearest-neighbor tij = 1.0, and 
lattice constant a = 1. 
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FIG. 1. A square lattice with nearest-neighbor hopping t 
and second neighbor hopping t' with the on-site Coulomb re- 
pulsion U and the nearest-neighbor repulsion V 



To determine the dominant gap function, we solve 
Eliashberg's equation with the FLEX approximation, 

WW = -§E r ( fc < fc ')G(fc')G(-fcW), ( 2 ) 



where <j> is the gap function, G Green's function, and T 
the pairing interaction with k = (k,w n ). The eigenvalue 
A, a measure of the pairing, becomes unity at T = Tc- 
For the calculation we take an N = 32 x 32 lattice, the 
temperature T = 0.03, and the Matsubara frequency 
for fcrmions -(2N C - l)nT < u n < (2N C - l)nT with 
N c = 1024. 

Esirgcn et al 31 ' 39 ' 40 have extended the FLEX method 
to general lattice Hamiltonians including the extended 
Hubbard model. Following them we introduce the pair- 
ing interaction, 

r B (fc,*')= E {§[^pX.pV^](fc-fc';Ar;Ar')e i(fc - Ar+fc '- A, " ) 

AT\AT" 

-i[FchXchKh] (k - fc'; Ar; Ar ' )e *(fc-Ar + fc'.Ar') 
+ iF s (0;Ar;Ary' fe ' Ar '- fe '' Ar > }, (3) 
for singlet pairing, and 

r t (fc,*')= E {-|^ P XspF sp ](fc-fc';Ar;Ar')e l ( fc ' Ar+fe '' Ar '» 

A7\Ar" 

-m*xM (k - fc'; Ar; Ar y(fc-Ar+fc'.Ar') 

+ iK(0;Ar;Ar> l ( fc ' Ar '- fc '' Ar )} (4) 

for triplet pairing. Here Ar(= 0,±x,±y) is null or 
nearest- neighbor vectors, 

Xs P = x/(i + K P x), 

Xch = x/(l + V ch x) 

are the spin and charge susceptibilities, respectively, 
where \ is the irreducible susceptibility, 

X(q; Ar; Ar') = e ik> < Ar - Ar ^G(k> + q)G(k>), (5) 

fe' 

and Vch(V sp ) is the coupling between density (magnetic) 
fluctuations, 

V ch (q;Ar;Ar') 

C U + 4[K cos(^) + V y cos(q y )], Ar = Ar' = 0, 
= <^ -14, Ar = Ar' = ±x (6) 

[ -V y , Ar = Ar' = ±y 

( -U, Ar = Ar' = 0, 
V ap (q; Ar; Ar') = < —V x , Ar = Ar' = ±x (7) 
1 -V v , Ar = Ar' = ±y, 

where q = (q, e n ) with e„ = 2nirT being the Matsub- 
ara frequencies for bosons, and V x = V y = V here. We 
have found that the q dependence of V^ p and V^h does 



not significantly affect T s and T t . Accordingly the peak 
position of Xch is almost the same as that for V^hXchVdi 
term in the expression for T. 

V s (0; Ar; Ar'), 14(0; Ar; Ar'), appearing in the last 
lines in eqs.(3,4) respectively, are constant terms involv- 
ing U and V, 



2U. 



V s (q;Ar;Ar') 



V t (q;Ar;Ar') 



V x e 



±iq u 



V, 



-V x e 

V v , 

-V y e ±iq y, 



Ar = Ar' = 0, 
Ar = Ar' = ±x 
Ar = -Ar' = ±x 
Ar = Ar' = ±y 
Ar = -Ar' = ±y 

V x , Ar = Ar' = ±x 
Ar = -Ar' = ±x 
Ar = Ar' = ±y 
Ar = -Ar' = ±y 



±iq. 



(8) 



(9) 



When the off-site interaction V is introduced all the 
vertices (V sp , V c h, V a , Vt) as well as the susceptibilities 
become (Z + 1) x (Z + 1) matrices for the lattice coordi- 
nation number Z{= 4 for the square lattice). 



III. RESULT 

A. A case study for Sr2Ru04 

We consider 7 band of Sr2Ru04 whose band filling is 
n = 4/3. To represent the shape of the Fermi surface ob- 
tained by ARPES 41 , we choose parameters n — 1.333 and 
next nearest hopping t' = 0.5. We first show the phase 
diagram against on-site Coulomb interaction U and the 
nearest neighbor Coulomb interaction V (Fig. 2). There, 
we have assumed that the dominant superconductivity 
pairing is the one that has the largest eigenvalue A in 
Eliashbcrg's equation calculated at T = 0.03. While the 
value of A for T = 0.03 is still much smaller than unity, it 
is difficult to extend the FLEX calculation to lower tem- 
peratures, given the complexity of the (extended Hub- 
bard) model. So this amounts to an assumption that the 
order in which A's appear for various pairing symmetries 
does not change for T — > 0. The charge density wave 
(CDW) is identified as the region in which the charge 
susceptibility (which has a peak at {tt,tt)) diverges at 
T = 0.03. 

In the result a triplet superconductivity phase (which 
does not exist for a single band when V is absent) is 
seen to appear just below the CDW phase for U > 9. 
The maximum eigenvalue of charge and spin suscepti- 
bilities and detailed behavior of A as a function of V 
for fixed U = 10 is shown in Fig. 3. From the figure 
we see that the triplet superconductivity becomes dom- 
inant when the charge susceptibility Xch is larger than 
the spin susceptibility Xsp- Figure 4 shows gap func- 
tions in fc-space, which shows that the triplet (which has 
the highest A) has a gap function oc sm(k x + k y )), while 
the two singlet gap functions have cos 2k x — cos 2k y and 
sin k r sin k v . A2 



We can trace back the reason why these gap functions 
are favored to the structure of spin and charge suscep- 
tibilities. Figure 5 shows the peak positions of spin and 
charge susceptibilities for V = 2.7. The peak positions 
for charge are around (±7r, ±7r), (±7r, ttt), while those for 
spin are ferromagnetic-like around (±j,±j), (±^, Tf )■ 
In Eliashberg's eq.(2)-(4), we see that for singlet the co- 
efficient of x S p (xch) term in T s are positive (negative). 
So the dominant gap function must change sign across 
the momentum transfer that points the peak position in 
the spin susceptibility, while the gap function must not 
change sign across the momentum transfer that points 
the peak position in the charge susceptibility. For triplet, 
on the other hand, the coefficients of x sp and Xch terms 
in r t are both negative. So the gap function should have 
the same sign across both the spin and charge peak po- 
sitions. 

When V is small, the spin susceptibility is much 
greater than the charge susceptibility, and the spin- 
fluctuation mediated pairing interaction for singlet is 
three times larger than that for triplet. This explains 
cos2fc x — cos2fcj,, which changes sign across (±^,±^), 
(ij, Tf )> becomes the dominant gap function. For large 
V, the charge susceptibility becomes dominant, and the 
pairing interaction becomes the same between singlet and 
triplet channels, i.e., T s and T t become similar in mag- 
nitude. Whether triplet is more favored than singlet de- 
pends on the other factors such as the number of nodes 
on the Fermi surface, which work unfavorably for pairing 
since the integration around the node in the right hand 
side of eqn.(2) cannot contribute to A. So a larger the 
number of nodes is basically unfavored. 

As another factor, we can find from Fig. 3 that the 
A's for the gap function sm(k x + k y ) and sinfcj, sinfcj, in- 
crease more sharply than that for cos 'Zkx cos2k y for 
V > 2, where Xch increases sharply. The reason should 
be that the first two have the same sign across the peak 
°f Xch around (±7r, ±7r) or (±7r, T 71 "); while the latter has 
slightly warped nodal lines. So there is a region exists 
where the gap function has the opposite signs across the 
peak of Xch- 

Returning to S^RuCU, our result in Fig. 3 indicates 
that if triplet can be dominant, the gap function should 
be sm(k x + k y ) which is degenerate with sin^ — k y ). 
The true gap function below T c should be a complex lin- 
ear combination, 



sm(k x + k y ) + ism(k x - ky), 



(10) 



which is more stable thermodynamically and breaks the 
time-reversal symmetry (Fig. 6). The absolute value of 
the gap function has minima on the Fermi surface along 
[100] and equivalent directions as depicted in Fig. 6 with 
open circles, which is consistent with the recent exper- 
iment of the specific heat in rotated magnetic fields. 13 
However, we also notice that the triplet region is very 
narrow, so identification of the pairing in Sr2Ru04 would 
require a precise determination of the parameters in that 



material. In addition, if we include the effect of a and (3 
band of Sr2Ru04 this triplet region may expand. 
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FIG. 2. Phase diagram against V and U with n — 1.333 
and t' = 0.5 for the 2D extended Hubbard model. 
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FIG. 4. The dominant gap function sin(fc : 
in k space for V = 2.7 (right panel), cos2fc a 
glet) for V — 2 (left), and sub-dominant singlet gap function 
sin fc^ sinky for V = 2.6 (center), each for U = 10.0. The ar- 
rows indicate the main scattering process mediated by charge 
fluctuation. 
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FIG. 3. Top: the maximum eigenvalue of Xch[k,0) charge 
(solid line) and Xsp(fcjO) spin (dotted) susceptibilities as a 
function of V for U = 10.0. Bottom: the eigenvalue A of 
Eliashberg's equation for triplet sin(k x + k y ) (solid line), sin- 
glet cos 2k x — cos 2k y (dotted) and sin k x + sin k y (dot-dash) 
as a function of V for U — 10.0. The CDW (gray) region is 
identified from the divergence in the charge susceptibility. 




FIG. 5. Color-coded plots of the spin susceptibility 
Xs P (fc,0) (left) and the charge susceptibility \ch(k, 0) (right) 
for U = 10.0, V = 2.7. 




FIG. 6. Contour plot of | sin(k x + k y ) +isin(k x — k y )\ along 
with the Fermi surface for T = 0.03, U = 10.0, V = 2.7. Open 
circles indicate minima of the gap function on the Fermi sur- 
face. 



B. Relation between the pairing symmetry and the 
Fermi surface 



Let us identify the relation between the pairing sym- 
metry and shape of the Fermi surface in the present con- 
text. For that purpose we can change the second neigh- 
bor hopping t' , which controls the warping (and even 



the topology) of the Fermi surface. In Fig. 7 we show 
A and the dominant pairing symmetry when t' is varied 
with fixed U = 10.0 and V = 2.7. Figure 8 depicts how 
the Fermi surface obtained by e£ + ReS(fe) = fi changes 
with if . The reason why those gap functions dominate 
can be understood in terms of the structure and value 
of the spin and charge susceptibilities. While the peaks 
in the spin susceptibility for V = correspond to the 
nesting vectors for the Fermi surface, the peaks for large 
V do not necessarily correspond to those. The result 
shows that the peaks change from (0.67r,7r), (7r,0.67r) for 
t' = - 0.3, (0.3tt, 0.3tt) for t 1 = 0.4 - 0.8, to (0, 0.5tt), 
(7r,0.57r) for = 0.9 - 1.0 as seen in Fig. 9. On the 
other hand the peak position in the charge susceptibility 
remains at (tt, tt) for the whole range of t' . The peak val- 
ues of Xch(fc, 0) and Xs-p(k, 0) are shown in Fig. 10, where 
the peak of Xch(k, 0) changes rapidly while x sp (fc, 0) has 
an almost constant peak value. 
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FIG. 7. The maximum eigenvalue A of Eliashberg's equa- 
tion for spin-triplet (solid line), spin-singlet (dotted) as a 
function of the second neighbor hopping t' for U = 10.0 and 
V — 2.7 with the dominant pairing symmetry indicated 



From the result we see that nearly antiferromagnctic 
spin fluctuations favor d x 2_ y 2 pairing for t' = — 0.3 
as in the high-T c cuprates, strong charge fluctuations fa- 
vor sin(k x + k y ) for t' = 0.3 - 0.6. For t' = 0.7 - 1.0 
where the spin susceptibility again exceeds the charge 
susceptibility, nearly ferromagnetic spin fluctuations fa- 
vor cos 2k x — cos 2k y which has many sign changes on the 
Fermi surface. 



Finally it is interesting to identify the reason 
why the charge fluctuation is sharply peaked around 
t' = 0.4. To elaborate the point, we first 
note the van-Hove singularities, which reside at 
(0,±tt), (±tt,0), (±&TCCOs(-t/(2t')),±a,rccos(~t/(2t'))) 
and (±arccos(-t/(2t')),Tarccos(-t/(2t'))) for the t-t' 
tight-binding model. So the charge susceptibility with a 
peak around (tt, tt) should be maximized when the Fermi 
surface approaches (0, ±7r) and (±7r,0), since the pair 
hopping with the momentum transfer (tt, tt) connects the 
van-Hove singularities, where the density of state is large. 
In the terminology of the previous work of ours 43 , the 
Fermi surface becomes "thick" around the van-Hove sin- 
gularities. 
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FIG. 8. Fermi surface (solid line) and nodes of the dom- 
inant gap function (dotted) for t' = 0.1 (top left), t' = 0.5 
(top right), t' = 0.7 (bottom left) and t' = 0.9 (bottom right) 
for U = 10.0 and V = 2.7. 



We have actually studied the relation between the 
charge susceptibility and the Fermi surface as shown in 
Fig. 11 to confirm that the charge susceptibility becomes 
maximum for t' = 0.41, at which the Fermi surface just 
touches the van- Hove singularities (0, ±7r), (±7r, 0) and 
the topology of the Fermi surface changes. We see in the 
figure that the Fermi surface becomes "thick" around the 
van-Hove singularities. The value t' = 0.4 is close to the 
t' for 7 band of Sr2Ru04, which may imply that triplet 
superconductivity induced by charge fluctuations may in 
fact be relevant there. 
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FIG. 9. Color-coded plots of the spin susceptibility 
Xs P (fc,0) for t' = 0.1 (left), t' = 0.7 (center), and t' = 0.9 
(right) for U = 10.0 and V = 2.7. 



IV. CONCLUSION 
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FIG. 10. The maximum eigenvalue of Xch(k, 0) (solid line) 
and Xs P (fc,0) (dotted) as a function of t' for U = 10.0 and 
V = 2.7. 
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FIG. 11. The t' dependence of the Fermi surface (solid 
lines), while dotted lines delineate E — Ef ±0.2 for U = 10.0, 
V = 2.7. 



We have studied the pairing symmetry in the one-band 
extended Hubbard model having the nearest neighbor 
Coulomb repulsion for square lattice. For the parame- 
ter which represents 7 band of Sr2Ru04, spin-triplet su- 
perconductivity whose symmetry is sin(/ca; ~\~ k y ) appears 
just below the CDW phase. So the true gap below T c 
is suggested to be a chiral s'm(k x + k y ) + isin(k x — k y ) 
which breaks the time-reversal symmetry. It would be 
interesting to experimentally detect the pairing state us- 
ing, for example, a theoretical prediction made in ref. 44. 
However, the region for the spin-triplet superconductiv- 
ity is very small, and the identification of the pairing in 
Sr 2 Ru04 will require a precise determination of parame- 
ters, and also a study of the effects of a and f3 bands. 

We have also examined the relation between the pair- 
ing symmetry and the shape of the Fermi surface by vary- 
ing t 1 . As t' increases, the shape of the Fermi surface 
changes drastically which changes the structure of x sp 
and x c h. This in turn changes the pairing symmetry 
as d x 2_ y 2 — * sa\{k x + k y ) — > cos2k x — cos2fc y . In addi- 
tion, Xch becomes especially large when the Fermi surface 
touches van-Hove singularities at which the Fermi surface 
changes the topology, which is achieved for t' = 0.41 near 
the parameter for 7 band of Sr 2 Ru04. 
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C. Real-space representation 

So far we have shown all the results in fc-space, but 
it is also heuristic to represent the pairing symmetry in 
real space as depicted in Fig. 12. From Fig. 7-10 we 
see that d x 2_ y 2 is suppressed as the (tt,tt) charge fluc- 
tuation increases, while from Fig. 3 sin(k x + k y ), d xy 
and cos2fc x COS 2fcy are favored as the charge fluctu- 
ation increases. This should be because strong (7r, 7r) 
charge fluctuations tend to arrange electrons diagonally 
on the lattice, and suppresses d x 2_ y 2 pairs across the 
nearest neighbors, but do not disturb the pairing such 
as sva.(k x + k y ), d xy and cos2fcz — cos2fcj, formed across 
more distant sites. 




FIG. 12. The pairing symmetry in real space. 
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